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Abstract 

Fukui et al. (2006) discovered two molecular loops in the Galactic center that are 
likely created by the magnetic flotation due to the Parker instability with an estimated 
fleld strength of ~150 fiG. Following the discovery, we present here a detailed study 
of the two loops based on NANTEN ^^CO(J=1-0) and ^^CO(J=1-0) datasets. 

The two loops are located in / = 355° - 359° and 6 = 0° - 2° at a velocity 
range of -20 - -180 km s~^. They have a projected total length of ~ 600 pc and 
heights of ~250 - 300 pc from the Galactic disk at a distance of 8.5 kpc. They have 
loop-like filamentary distributions of 30 pc width and show bright foot points in the 
^^CO emission in the edges of the loops at 6 ~ 0.8° - 1.0°. These foot points are 
characterized by velocity dispersions of 50 - 100 km s~^, much larger than those in 
the Galactic disk, supporting that the loops are located in the Galactic center within 
~1 kpc of Sgr A*. The loops also show large-scale velocity gradients of ~30 - 50 km 

per ~100 pc. 

We present an attempt to determine geometrical parameters and velocities of the 
loops by assuming that the loops having the same size are expanding and rotating at 
a constant radius R. The analysis yields that the loops are rotating at 50 km s~^ and 
expanding at 150 km s~^ at a radius of 670 pc from the center. 
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We find some new features on smaller scales in the loops; first, the tops of loops 
1 and 2 exhibit helical distributions in certain velocity ranges and, second, one of 
the foot point of loop 2 shows a sub-structure which mimics a smaller loop. These 
features possibly offer further indications for the strong magnetic field in the loops 
because of the similarity with the solar features resulting from magneto- hydrodynamic 
instabilities. 

The loops have counterparts of the H I gas indicating that atomic gas also consists 
the magnetic loops. The H I gas with weaker intensity fills the inner part of the loops 
where there is no sig nificant ^^CO and ^^CO emission. The ^^CO, ^^CO and Hi 
datasets were used to estimate the total mass and kinetic energy of loops 1 and 2 
to be ~ 5.5(±2.4) x IO^Mq and ~10^^ ergs, respectively. We also examined the far 
infrared distributions of IRAS toward the loops and found that the dust emission is 
well correlated with the gas distribution. The atomic hydrogen mass estimated from 
the dust emission is consistent with the above estimates. 

Based on the new physical parameters, we present a reinforced argument that 
the huge size, velocity dispersions and that the kinetic energy are consistent with 
the magnetic floatation origin due to the Parker instability but is difficult to be 
explained by multiple stellar explosions. We also present detailed comparisons with 
the numerical simulations of magnetized nuclear disk by Machida et al. (2009) and 
Takahashi et al. (2009) and note that the model predictions show excellent agreements 
with the observations. 

Key words: Radio lines: ISM — ISM: clouds — ISM: magnetic fields — magnetic 
loops 

1. Introduction 

Molecular clouds are the sites of star formation and the distribution and dynamics of 
molecular gas should be crucial in the evolution of galaxies. The evolution of molecular gas 
may be significantly different in the central region of a galaxy from that in the disk because the 
stellar gravitational field is considerably stronger in the center due to increased stellar density. 
The strong stellar gravitational field makes the pressure much higher in the central region than 
in the disk and should influence considerably the gas dynamics. It is important to understand 
the physical conditions of molecular gas and their astronomical consequences in the Galactic 
center in our continuing efforts to elucidate galactic evolution. 

Previous studies of the Galactic center show most of the molecular gas in the Galactic 
center is concentrated in the inner 300 pc which is called the gCentral Molecular Zoneh (here- 
after CMZ, e.g., Morris & Serabyn 1996). CMZ including Sgr A and Sgr B2 molecular clouds 
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has been a region intensively studied in the Galactic center (e.g., Fukui et al. 1977; Glisten 
& Henkel 1983; Scoville & Solomon 1974). There are also molecular features with weaker 
intensities outside the CMZ up to nearly 1 kpc from the center. They include clump 1 ~ 
355°, 0°), Clump 2 ~ 3°, 0°) and 5-deg feature {l,b ~ 5°, 0°) (Bania 1977) and at least 
several weak CO emission features detected at a lower resolution of 8'. 8 (Dame et al. 1987; 
Bitran et al. 1997) although these features outside the CMZ received little attention so far. 

Fukui et al. (2006; hereafter Paper I) discovered that two molecular loops are located 
toward 355.5° < / < 358° and 0° < 6 < 2° by analyzing the NANTEN Galactic plane survey 
dataset of the ^^00(7=1-0) emission. These loops have heights of ~ 2 degrees from the 
galactic plane and the projected lengths of 3 - 4 degrees. The loops have two foot points that 
are brightest spots in the CO emission on each end at Galactic latitudes around 0.8 degrees. 
The velocity dispersions of the molecular gas in the loops are ~50 - 100 km s~^ at the foot 
point of loops, which are characteristic to the molecular gas in the Galactic center and are 
much larger than that of the disk molecular clouds whose velocity widths are less than ~ 10 
km s~^.@Most recently, Fujishita et al. (2009) have discovered loop 3 which is located in the 
same direction with loops 1 and 2 at a positive velocity range of Vlsr = 20 - 200 km s~^. Loop 
3 shows a loop-like outer shape and similar height and length with loops 1 and 2 and it has 
two foot points with broad velocity dispersions. 

These three loops are interpreted by magnetic flotation caused by the Parker instability 
and Paper I presented the first observational evidence for the Parker's model in the galactic 
scale. The differential rotation in the nuclear disk creates toroidal magnetic field where the 
molecular gas is frozen-in to the magnetic field; the molecular gas in the disk is ionized at an 
ionization degree of 10^*^ - 10^^ by cosmic ray protons and the rate is high enough for the 
frozen-in condition that are described by the ideal magneto-hydrodynamics. The molecular gas 
layer in the disk is dynamically supported by the toroidal field against the stellar gravity in the 
z direction. This configuration is unstable to disturbances and the gas rises up as a loop by the 
Parker instability (Parker 1966). The gas inside the loop then flows down to the disk by the 
stellar gravity and the down flowing gas forms shock fronts above the disk (Matsumoto et al. 
1988). This shock interaction forms large velocity dispersions at the foot points of loops, heating 
and compressing the gas by the shock fronts (Paper I). In particular, the most recent global 
numerical simulations that deal with the 2kpc-radius nuclear disk by Machida et al. (2009) 
has revealed general properties of the magnetic floatation loops over the entire disk including 
that the one-armed non-axisymmetric (m = 1) mode tends to dominate the gas distribution. 
Subsequently, Takahashi et al. (2009) carried out two-dimensional simulations of a local part 
of the Galactic center disk. They assumed a stratified disk consisting of cool layer (T ~ 10'^ K), 
warm layer (T ~ 10^ K), and hot layer (T ~ 10^ K) to confirm the observed galactic loops, and 
found that the flat-top loop efficiently accumulate the gas to make a dense layer at the top of 
the loop. 
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In the CMZ, the molecular gas shows high temperatures and violent motions. The 
temperature is derived to be ~30 - 300 K using CO, NH3, H2, H3' lines (Martin et al. 2004; 
Hiittemeister et al. 1993; Rodriguez- Fernandez et al. 2001; Oka et al. 2005; Nagai et al. 2007). 
The velocity dispersion is 15 - 50 km s~^, which is significantly larger than that in the molecular 
clouds in the Galactic disk (e.g., Morris & Serabyn 1996; Giisten & Philipp 2004). The causes 
of the high temperature and large velocity dispersion have been the longstanding puzzles in the 
last few decades. The magnetic field is as strong as 0.1 - 1 mG as indicated by Radio Arc and 
other non-thermal filaments which are distributed within 1 degree of the center, often vertically 
to the Galactic plane (Yusef-Zadeh et al. 1984; 2004), and perhaps by the infrared double helix 
(Morris, Uchida & Do 2006). The magnetic fiotation model has a potential to offer a coherent 
explanation of the origin of the high temperature and large velocity dispersion of molecular 
gas in the Galactic center as argued in Paper I because the high velocity motions are natural 
consequences of the gas motion along the magnetic loops which form shock fronts around them 
most significantly at the foot points. 

Following the discovery of loops 1 and 2, it is important to clarify details of the loops 
including additional indications of the shock interaction. We here present a detailed analysis of 
loopsl and 2 by using the NANTEN Galactic plane survey, the GPS ^^CO and ^^CO datasets. 
Section 2 describes the ^^CO and ^^CO(J=1-0) datasets. Detailed observational properties of 
the magnetic fiotation loops are presented in section 3. Section 4 offers discussion and section 
5 summarizes the paper. 

2. NANTEN Datesets 

Basic parameters of the observations are summarized in Table 1. 
2.1. ^^CO 

^^CO(J=1-0) dataset toward the Galactic center region was taken with the NANTEN 
4-m radio telescope at Las Campanas Observatory in Chile during the period from March 1999 
to September 2001. The half-power beamwidth (HPBW) was 2'. 6 at 115 GHz, the frequency 
of ^^CO(J=1-0). The front end was a 4 K cryogenically cooled Nb superconductor-insulator- 
superconductor (SIS) mixer receiver (Ogawa et al. 1990) that provided a typical system tem- 
perature of ~280 K in the single-side band, including the atmosphere toward the zenith. The 
spectrometer was an acousto-optical spectrometer (AOS) with 2048 channels. The frequency 
coverage and resolution were 250 MHz and 250 kHz, corresponding to a velocity coverage of 
650 km s~^ and a velocity resolution of 0.65 km s~^, respectively, at 115 GHz. The intensity 
calibration was made by using the ambient temperature load. The absolute antenna tem- 
perature was calibrated by observing rho-Oph East [R.A.(1950) = 16^29™20f9, Dec. (1950) = 
— 24°22'13"] every 2 hours, whose absolute temperature was assumed to be 15 K. The telescope 
pointing was measured to be accurate to within 20" by radio observations of planets in addition 



4 



to optical observations of stars with a CCD camera attached to the telescope. The observed 
region was 240 square degrees toward -12° < / < 12° and -5° < 6 < 5° at a grid spacing of 4', 
corresponding to 10 pc at a distance of the Galactic center, 8.5 kpc. In total, 54,000 positions 
were observed. The integration time per point was 4 - 9 s, resulting in typical r.m.s. noise 
fluctuations of 0.36 K at a velocity resolution of 0.65 km s~^. 

2.2. ^^CO 

^^CO(J=1-0) dataset toward the Galactic center region was taken with the same instru- 
ment as the ^^CO(J=1-0) observation in October 2003. The half-power beam width (HPBW) 
was 2'. 7 at 110 GHz, the frequency of ^^CO(J=1-0). The same receiver provided a typical 
system temperature of ~100 K in the single-side band, including the atmosphere toward the 
zenith. The same spectrometer provided the frequency coverage and resolution, 250 MHz and 
250kHz, corresponding to a velocity coverage of 620 km s~^ and a velocity resolution of 0.62 km 
s~^, respectively, at 110 GHz. The absolute antenna temperature was calibrated by observing 
rho-Oph East [R.A.(1950) = 16''29'^20f9, Dec.(1950) = -24°22'13"] every 2 hours. Absolute 
radiation temperature of rho-Oph East was assumed to be 10 K. The observed region was 28 
square degree toward -6° < / < 8° and -1° < 6 < 1° at a grid spacing of 2'. In total, 25,200 
positions were observed. The integration time per point was 10 - 15 s, resulting in typical 
r.m.s. noise fluctuations of 0.20 K at a velocity resolution of 0.62 km s~^. 

3. Results 

3.1. CO Distribution of Loops 1 and 2 

Figure 1 shows the large-scale distribution of ^^CO intensity in an area of / = 350° - 10° 
and b = -5° - 5° in Vlsr= -300 - 300 km s~^. Figure la shows the total integrated intensity 
distribution of ^^CO. In Fi gure la most of the emission along the galactic plane is the 3-kpc 
expanding arm and the CMZ is dominant in / = 359° - 2°. Locations of loops 1 and 2 and loop 
3 toward / = 355° - 359° are shown by boxes with broken lines in Figure la and Figure lb. 
In Figure lb, the longitude velocity diagram, the molecular gas in the Galactic center region 
is seen as broad features whose velocity extents are a few tens km s~^ or more. The narrow 
features near km s~^ are clouds outside the Galactic center and the 3-kpc arm is seen as a 
narrow feature at ~-50 km s~^ at / = 0° with a velocity gradient of ~8 km s~^/deg. 

Figure 2 shows the integrated intensity distributions for the negative and positive velocity 
ranges. Figure 2a and 2b show the ^^CO distributions in the negative velocity range in the l-b 
plane and b-v plane, respectively, and loops 1 and 2 are indicated in Figure 2a. Figure 2c and 
2d show the ^^CO distributions in the positive velocity range in the l-b plane and b-v plane, 
respectively, and loop 3 is indicated in Figure 2c. 

Figure 3 shows close-up images of ^^CO and ^^CO for loops 1 and 2 for two velocity 
ranges, from -180 to -90 km s~^ and from -90 to -40 km s~^, in an area of / = 355° - 359° 
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and 6 = 0° - 2° shown by a rectangle in Figure 2a. Figures 3a-d are for loop 1; Figures 3a 
and 3b are the l-b diagrams of ^^CO and '^^CO, respectively, and Figures 3c and 3d are the l-v 
diagrams of ^^CO and ^^CO, respectively. Fi gures 3e-h are for loop 2; Figures 3e and 3f are 
the l-b diagrams of ^^CO and ^^CO, respectively, and Figures 3g and 3h are the l-v diagrams 
of ^^CO and ^^CO, respectively. 

Figure 3a shows that loop 1 is a remarkable filamentary feature of ~30 pc width and 
~300 pc projected length, being elevated from the plane by ~200 pc. These parameters are 
all assumed the distance to the Galactic center as 8.5 kpc. It shows a bright spot at {l,b) ~ 
(356.0°, 1.0°), which has strong intensity gradients toward the plane and also toward the east. 
We also recognize a less enhanced peak toward {l,b) ~ (357.4°, 0.8°). We shall call these two 
peaks gfoot pointsh following Paper I. Figure 3c shows that these two foot points have enhanced 
extents of ~50 km s"^ or more. The ^'^CO emission is significant only toward the foot point at 
(356.0°, 1.0°) while the observations does not cover the region above b = 1°. 

Figure 3e shows that loop 2 is an outstanding feature with a width of ~30 pc and a 
length of ~400 pc and is elevated by ~300 pc from the plane. A feature at b ~ 1.4° extending 
toward the center is part of loop 1. Loop2 shows two bright spots with enhanced extents of 
~60 km s-i toward {l,b) ~ (355.4°, 0.8°) and (356.1°, 0.8°). The ^^CO emission is detected 
toward these two spots (Figure 3c, d) although the coverage is limited to b less than 1°. We 
shall call these two spots gfoot pointsh of loop 2 (Paper I). 

Figure 4 shows typical ^^CO and ^^CO line profiles toward loops 1 and 2. These show 
the line profiles are not Gaussian but have shapes of peaked triangles. We note that the profiles 
toward the foot points show asymmetry with a strong skew. The velocity span is as large as 40 
km s~^ towards the loop top and becomes even larger as 100 km s~^ towards the foot points. 

Figures 5a and 5b show distributions of the averaged velocity in the loops where the 
averaging was made in a 5 arcmin Gaussian beam. The larger beam size was taken in order to 
suppress small scale irregularities. We roughly estimate velocity gradients of ~0.2 km s^^ pc~^ 
and ~0.35 km s"^ pc~^ along the loops for loops 1 and 2, respectively. 

In order to show more details of the velocity distributions we show a series of velocity 
channel distributions every 10 km s~^ as 12 panels in Figure 6. We identify the foot points of 
loops 1 and 2 towards / = 355°, 357° and 359° b = 0.8°. From -120 to -90 km s"^ we find a 
vertical feature toward / = 356° at b from 1.4° to 1.6°. This feature between loops 1 and 2 does 
not seem to be part of the loops. 

In addition, we note three regions show peculiar distributions that are not part of simple 
loops as shown in Figures 7, 8 and 9. 

Figure 7 shows the details of the top of loop 1 every 4 km s~^ from -129 to -107 km 
s~^. The panels from -121 to -113 km s~^ show S-shaped distributions in a range / = 357.0° - 
357.4° and those from -115 to -113 km s"^ also show a similar feature in a range from (356.8° 
- 357.3°, 0.9°). 
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Figure 8 shows the top of loop 2 every 4 km s~^ from -86 to -66 km s~^. A pecuhar 
peaked shape is found from / = 355.4° to 355.8°. and b = 1.8° to 2.1°. Figure 9 shows the area 
in the north of the foot point. The panels from -50 to -42 km s~^ show a vertical feature at 
/ = 355.6° and the panels from -44 to -36 km s~^ indicate a bridge-like feature that connects 
this vertical feature and the local intensity peak at (/,6) = (355.0° - 355.2°) in panel from -44 to 
-40 km s~^. In addition, the panel from -34 to -30 km s~^ shows another vertical feature at I 
= 355.2° and b = 1.4° to 1.8°. 

3.2. Comparison with H/ 

Figure 10 shows a comparison between ^^CO and Hi where the Hi data were taken with 
the Parkes 64m telescope having 16 arcmin beam (McClure-Griffiths et al. 2005). Figures 10a 
and 10b are for loop 1 and Figures 10c and lOd for loop 2 in the l-b and l-v diagrams. The 
Hi distribution seems similar to ^^CO in the both loops indicating that the loops accompany 
atomic components also, while the Hi seems enveloping the CO. It is notable that the inside of 
loop 1 shows weak Hi emission (Figure 10a). Loop 2 is not so clearly recognized as loop 1 but 
the western half of loop 2 shows a possible Hi counterpart in Figure 10b. The outside of loops 
1 and 2 shows little sign of Hi emission. Figure 11 shows b-v cuts of ^^CO and Hi. Loop 1 is 
clearly seen from v ~ -150 to -100 km s~^ around b ~ 1.0 (Figure 11a) with the Hi associated 
with the CO (Figure lib). 

3.3. Comparison with Dust emission 

In order to test if the loops in CO and H I are seen in the dust emission we have compared 
dust emission obtained with IRAS at 60 and 100 /im, where we did not use shorter wavelength 
data because of the possible contamination by the zodiacal light. Figures 12a and 12b show 
comparison between the total integrated intensity of ^^CO from -300 to 300 km s~^, including 
loops 1, 2 and 3, and IRAS 60 /xm and 100 fim emissions, respectively. We find that the 
enhanced CO emission above 6 = 1° in a range of / = 355° - 358° shows a good coincidence 
with the dust emissions; the hole of the CO emission at / = 356.5° and the vertical feature at 
/ = 355.3°, the western part of loop 2, are recognized well particularly at 100 fim. Figures 13a 
and 13b, a similar comparison with the Hi convolved to the same beam size, show that the dust 
emissions show similar elevated distribution as in Figure 12. Figure 14 shows four panels of 
CO in different velocity intervals and suggests that both loops 1 and 3 contribute significantly 
to the main dust emission at b above 1° and / = 355° - 358°, in addition to the western half of 
loop 2 at / ~ 355.3°. 

We estimated the total hydrogen column density by combing the Hi and CO emissions 
(Figure 15a) and from the dust emission (Figure 15b). Here we describe the total column 
densities estimated by gas and dust as A^(Hl)gas and A^(Hl)dust, respectively. We used the 
following relationship to convert the dust emission into A^(Hl)dust by assuming the gas to dust 
mass ratio, Rgd, of 100; 
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iV(Hl).„.. = i (^) (cm-) (1) 

where a is the grain radius, p is the grain density in g cm~^, Qioo is the grain emissivity at 
100 fim, and mn is the atomic hydrogen mass. We use an average value for (ap/Qioo) of 3.2 
X (1000/A^m)~^5 assuming /3 = 2, at wavelength A for a mixture of graphite and sihcate grains 
(Hildebrand 1983; Agladze et al. 1996). The optical depth at wavelength A, tx, is represented 
as following; 

fx 

-B;^(A,Tdust) 

where fx is flux density at A, and B^{X,Tdnst) is Planck function. We estimated dust tempera- 
ture, Tdust; by comparing flux densities of 60 fim and 100 fim at each observed point shown in 
Figure 15 with following equation; 



/eo _ /60\-Y^60 

/loo viooy [q 



100 , 



exp {he /XioohTdnst) - 1 



60 



n 



100 



(3) 



exp{hc/ XeohT^^st) - 1 

where h is Planck constant, and c is the light speed, and kb is Boltzmann constant, and Qqo 
and f2ioo are the solid angles at 60 /xm and lOO/im, respectively. We assume that figo — ^loo- 
The relationship between A^(H l)gas and both the CO and H I emissions is given as follows; 

iV(Hl)gas = (l^(CO) X 1.6 X 10=^° X 2) + (l^(Hl) x 1.823 x 10^^) (cm^^) (4) 

where W is the integrated intensity in K km s~^ units. We use the conversion factor from 
I^(CO) to the molecular hydrogen column density N{B.2), what we call "X-factor", of 1.6 x lO^'^ 
cm~^/(K km s~^). We shall derive this value in the next subsection. 

Figure 16 shows the scatter plot of the two values of A^(H l) for the region shown in Figure 
15. The iV(Hl) values show a good correlation with a correlation coefficient of logio(iV(H)dust) 
= 0.76 X logio(A^(H)gas) + 5.47 and support the physical association between the dust emission 
and the gaseous loops. 

3.4- Mass estimate 

We made mass estimates of loops 1 and 2 by using three methods; i) molecular mass by 
using the X-factor newly derived from the present ^'^CO and ^^CO ratio, ii) Hi mass, and iii) 
dust mass. 

Total molecular mass is estimated using the X-factor, which is an empirical conversion 
factor from ^^CO integrated intensity to molecular hydrogen column density: 

X = N{R2)/WC^C0) (cm-V(Kkms-i)) (5) 

Here we use the integration ranges same those in Figure 3; i.e., -180 to -90 km s~^ for loop 1 
and -90 to -40 km s~^ for loop 2. The mass of the molecular gas is calculated using following 
equation: 
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MC^CO) = umuY.[D^^N{ii2)] (6) 

where, u is the mean molecular weight which is assumed to be 2.8, D is the distance of loops 1 
and 2, 8.5 kpc, and Q is the solid angle subtended by a unit grid spacing 4' x 4'. Two values 
of the X-factor, 0.24 x lO^o cm-V(K km s"!) (Oka et al. 1998) and 0.6 x lO^o cm'^ K km 
s~^ cm~^/(K km s~^) (Onishi et al. 2004) were previously derived in the Galactic center by 
using the virial mass. The dynamical state is however not in dynamical equilibrium when the 
magnetic flotation is operating and may be considerably different from the virial equilibrium. 
We therefore derive newly an X-factor by using the LTE mass derived from the present ^^CO 
and ^^CO data (Figure 3). 

Although ^^CO observations are confined to 6< 1°, ^^CO column densities are calculated 
by assuming the local thermodynamic equilibrium (LTE) condition to estimate the lower limit of 
molecular mass. The optical depth of ^^CO, r(^'^CO), was calculated using following equation: 



r(i3co) = -ln 



(7) 



5.29 X (J(Tex)- 0.164) 

where, T^(^^CO) and T^^ are the radiation temperature and the excitation temperature of 
^^CO, respectively. J(T) is defined as J(T) = 1/ [exp(5.29/T) - 1]. NC^CO) was estimated 
from: 

where, AV^ is the ^^CO linewidth. A^(^'^CO) is calculated for each observed position which 
matches the observed position in ^^CO. In this study. We assume T^x of 40 K, which is typical 
kinetic temperature estimated by a LVG analysis with multiple CO excitation lines in the foot 
point found / ~ 356.0° reported by Torii et al (2009b). Then, we adopt the ratio [H2]/[^^C0] = 
10"^ (Lis & Goldsmith 1989) to convert iV(^^CO) into iV(H2). The results are shown in Figure 
17 which show that X-factor is estimated as 1.6(±0.7) x 10^'^ [cm~^/(K km s~^)]. This value 
is close to those derived from gamma ray observations for the inner disk 1.6 X 10^° (Hunter et 
al. 1997), 1.9 X 10^° (Strong & Mattox 1996) and is an order of magnitude larger than those 
derived from the virial theorem. By using this X-factor, the total molecular masses of loop 1 
and 2 are calculated as 2.1(±0.9) x lO^M©, and 3.4(±1.5) x IO^Mq, respectively. 

Hi column densities are calculated assuming that the Hi 21 cm line is optically thin 
and Tg hu/k. The following equation is used: 

iV(Hl) = 1.82 X 10^^ X W{m) = 1.82 x 10^^ y T^dv (cm'^) (9) 

The masses of the atomic gas of loops 1 and 2 are then derived as 2.4 x 10^ Mq, including the 
mass of 0.7 xIO^Mq inside of the loop, and 2.5 x 10^ Mq, respectively. 

We also find that the dust emission likely corresponds to the loops at 60 /im and 100 fim 
of the IRAS data as shown in Figure 12. It is not certain what fraction of the IRAS emission 
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is from loops 1 and 2 because there is another loop, loop 3, in the same direction. We find 
above that the total dust mass toward loops 1, 2 and 3 is well proportional to the dust mass 
as expressed by a linear relationship (Figure 16) between the far infrared emission and the CO 
and H I emissions, and the dust emission is well explained as associated with the molecular and 
atomic gas in the loops for the nominal gas to dust mass ratio of 100. Here we assume that the 
dust emission is proportional to WC^CO). WC^CO) of loops 1 and 2 are 2.4 x 10^ K km s ^ 
and 3.9 x 10^ K km s~\ respectively, and W^^CO) of loop 3 is 1.4 x 10^ K km s"^ The relation 
ship in Figure 16 then indicates the total dust masses of loops 1 and 2 is 2.1 (±0.9) x 
and 3.4(±1.5) x IO'^Mq, respectively. 

The all estimated masses for loops 1 and 2 are summarized in Table 2. The total masses 
of loops 1 and 2 are reached to 2.3(±0.9) x lO^M© and 3.7(±1.5) x IO^Mq, respectively, and if 
we assume the uniform velocity dispersion of 30 km s~^ in the loops, the total kinetic energy 
of each loop amounts to ~2-3 xlO^^ erg. 

4. Discussion 

4.I. Magnetic flotation 

Parker (1966) presented the pioneering work on the Parker instability in order to explain 
cloud formation of the foot point of magnetic loop in the galaxy. Matsumoto et al. (1988) and 
Horiuchi et al. (1988) made detailed studies of the Parker instability analytically with a linear 
analysis and numerically including nonlinearity. Their results indicate that the fundamental 
parameters are the Alfven speed Va = -B/v^vrp and the pressure scale height H. The height 
and length of a loop is given as a few times H and several times H, respectively. A typical 
timescale is given by the ratio H/Va- The molecular and atomic gas is generally frozen-in 
to the magnetic field lines even with weak ionization at ionization degrees of 10~® - 10~^ 
and the magneto-hydrodynamical (MHD) description holds well in the interstellar space. An 
outstanding observational property of a magnetic fiotation loop is that the falling down gas 
becomes often supersonic and forms shock fronts at the both ends of the loop, where gas 
density and velocity dispersion become enhanced, and these compressed regions appear as foot 
points of enhanced density (Matsumoto et al. 1988). 

Two molecular loops, loops 1 and 2, were discovered in the Galactic center and an 
interpretation was presented that loops 1 and 2 are created by the Parker instability in the 
magnetized nuclear disk with a field strength of ~150 /xG (Paper I). These authors estimated 
the field strength by assuming energy equi-partition between the turbulent gas motion and the 
magnetic field. They also made two-dimensional numerical simulations of loops 1 and 2 and 
showed that the two loops are successfully reproduced by the magnetic flotation mechanism. 
The gas inside the loop flows down to the disk under the influence of the disk gravity, forming 
a uniform velocity gradient along the loop. This down flowing gas collides with the disk and 
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often forms shock fronts with enhanced turbulent motions. The two loops in fact show foot 
points on their both ends as is consistent with the theoretical prediction. It is likely that the 
gas in the foot point is heated and compressed by the shock fronts. Therefore, the magnetic 
flotation has a potential to offer a coherent explanation both on the large velocity dispersions 
and high temperatures of the nuclear gas disk. It is however to be noted that these simulations 
were local in the disk, not including the effects of rotation. 

Subsequently, Machida et al. (2009) carried out three-dimensional global numerical sim- 
ulations of the nuclear gas disk and showed that the magneto-rotational instability coupled with 
the Parker instability works to create more than a few 100 loops over a Ikpc-radius nuclear 
disk, where the axially symmetric Miayamoto-Nagai potential (Miyamoto & Nagai 1975) was 
adopted as the stellar gravitational field. Machida et al. (2009) have shown that one-armed 
non-axisymmetric density pattern of m = 1 mode is developed in the nuclear disk. In the 
half of the azimuthal area of the nuclear disk having lower density, magnetic pressure tends to 
become stronger compared with gas pressure and prominent magnetic loops tend to be formed 
preferentially on this side of the disk rather than in the other half having higher density. So, 
the predicted asymmetry of the global distribution of the loops in fact seems to be consistent 
with the observations that about three- fourths of the dense molecular gas, CMZ, is distributed 
in the positive Galactic longitudes and that three loops are distributed in the negative Galactic 
longitudes. 

4-2. Alternative explanations 

Supershells created by multiple massive stellar explosions may be a possible explanation 
alternative to the magnetic floatation. In the Galaxy supershells associated with molecular 
clouds have been studied in about ten regions; they includes the Carina Flare (Fukui et al. 
1999; Dawson et al. 2008), Gum nebula (Yamaguchi et al. 1999) and other molecular shells 
(Matsunaga et al. 2001). The typical sizes of these shells are from 100 to 500 pc and the total 
velocity dispersions are 20 km s~^ corresponding to the total kinetic energy of ~ 10^^ ergs. As 
derived in the previous section, the total kinetic energy of the loops is 10 times larger than this 
energy involved in molecular supershells and, more importantly, the velocity spans of the loops 
are from 50 to 100 km s~^ in total, considerably larger than that of supershells. Considering 
the high pressure in the Galactic center region, we infer that supershells in the Galactic center 
may have even smaller velocity dispersions. In addition, we see no indications of the stellar 
remnants in the central regions of the loops that could have driven the shells. We therefore 
conclude that the loops are not explained by supershells. 

We shall here discuss the connection between the bar like gravitational potential and the 
magnetic flotation. The bar like gravitational potential is considered to be a viable mechanims 
to explain the radial motions of the gas in the central few 100 pc of the Galaxy (Binney et 
al. 1991). Numerical simulations were used to find that the flow of gas in the Galactic center 
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is dominated by a bar that has corotation at r = 2.4 ± 0.5 kpc, which is viewed at an angle 
of 16 ±2° from its major axis (Binney et al. 1991). From the structure of the Hi terminal 
velocity envelope, it is deduced that the central mass density scales as p oc r"^-"^^ out to at 
least about 1.2 kpc along the bar's major axis. Consequently, the circular velocity curve is 
rising significantly through the radius range where a naive analysis of the tangent velocity 
leads to a falling rotation curve. Obviously, the bar like potential is relevant to create large 
scale non-radial motions but is not necessarily exclusive in the magnetic flotation picture. In 
order to incorporate the effects of the bar like potential, it is required to develop a model which 
adopts the bar potential instead of the Miyamoto-Nagai potential in the global simulation of 
the magnetized gas disk. It is naturally expected that magnetic flotation loops are also created 
in such calculations because the basic physics remains the same in connection with the Parker 
instability. 

4.3. Model fitting 

The numerical simulations by Machida et al. (2009) show that the magnetic loops tend 
to be distributed an approximately constant radius from the center. Fujishita et al. (2009) 
derived the geometrical and kinematic properties of the newly discovered loop 3 by assuming 
the constant radius. We shall here estimate the geometrical and kinematical parameters such 
as radius, rotation and radial velocities. 

First, we estimate the geometrical parameters of the loops. For simplicity we assume 
that the two loops have the same radius and physical length and shape and subtend the same 
angle with respect to the center. By considering triangles the corners of which are the Galactic 
center, the Sun, and a foot point of loops (the left foot point of loop 1: Iq ~ 357.5°, the right 
foot point of loop 1 and the left foot point of loop 2: /i ~ 356.0°, and light foot point of loop 
2: I2 ~ 355.5°), we can derive simultaneous equations by applying sine theorem to the three 
triangles; 

R = ^^rPh (10) 

sm(^o + ^o) 
i?osin(/i) 

sin(^o + ^ + /i) ^ ^ 

_Rosin(^2)_ .^21 
sin(^o + 2^ + /2) ^ ' 

where R is the radius, and Rq is the distance to the Galactic center, 8.5 kpc, and 6 is the 
angles subtended by the arcs between the foot points at the Galactic center, and ^0 is offset. 
Schematic view of these parameters are shown in Figure 18. We then solve the equations and 
derive R, 6q, and 6 as ~670 pc, ~27.5°, and ~31°, respectively. Therefore, the location of 
the left end of loop 1 is at 6q ~ 27.5°, the right end of loop 1 and the left end of loop2 is at 
6*0 + ^ ~ 58.5°, and the right end of loop 2 is at ^0 + 26* ~ 89.5°. The projected length to the 
disk of each loops, L, is 360 pc, which is consistent with the results of Takahashi et al. (2009). 
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Next, we try to estimate the kinematical parameters of the loops. It is uncertain if the 
loops have radial motions. Some other features like the expanding molecular ring (EMR) are 
supposed to be expanding (Sawada et al. 2001; Morris & Serabyn 1996). We shall for simplicity 
assume that the loops are rotating and expanding from the center although the expansion is 
not yet confirmed in loops 1 and 2. We test the same method as discussed by Fujishita et al. 
(2009). All the molecular gas of loop 1 and 2 has negative velocities. The velocity with respect 
to the LSR, Vlsr, is then expressed as follows; 



where, Vj-ot is the rotational velocity of the disk, V^xp is the expansion velocity, and Vgun is the 
rotation velocity of the LSR about the center of 200 km s~^. We estimate that V^xp is 141 
km s~^ and V^t is 47 km s~^. These parameters are listed in Table 3 for the three loops, and 
Figure 18 shows schematic face-on view of these loops. The radius of loops 1 and 2 is somewhat 
smaller than that of loop 3. 

4-4- Foot points; evidence for shock formation 

Foot points are characteristic features in the magnetic flotation loops. The general 
observational properties of the foot points are summarized as follows while the eastern foot 
point of loopl seems to be under-developed compared with the other three; 1) The size and 
mass of the gas in the foot points are ~20 pc and ~10^ Mq. 2) The spatial distribution shows 
a sharp intensity drop toward the plane 3) Velocity dispersions are as large as 50 - 60 km s~^. 
4) The line profiles are of triangle shape with some skewness. 

The two foot points in loop 2 seem quite consistent with the model prediction by 
Matsumoto et al. (1988) because they are located nearly at the same altitude of ~ 150 pc 
above the galactic plane. The foot points also show strong intensity gradient toward the plane, 
suggesting shock fronts. The numerical simulations by Machida et al. (2009) show that the 
large velocity dispersion is excited by the energy released by the falling-down gas. A LVG 
analysis of the CO J=l-0, 3-2, and 4-3 transitions indicates that the typical gas temperature 
is about 50 K in the foot points of loops 1 (Torii et al. 2009b). We also note that the high 
density in the foot points, ~ 10^ cm~'^ may offer a site of future star formation. 

4-. 5. Formation of molecular loops 

Generally speaking, it is easier for lower density gas to be lifted up rather than dense 
molecular gas by the magnetic flotation and it may seem odd that the dense molecular gas 
forms the loops. Observations suggest that such a loop may be initially formed within the disk 
as an elongated molecular gas tube and then rises up as observed, since the eastern part of loop 
1 below the foot point may indicate such 

Paper I made numerical simulations of loops for two-temperature, hot and cool, layers 




(13) 



13 



and showed that magnetic flotation works to form loops as observed. In order to reproduce the 
observed cool molecular loops in the Galactic center; Takahashi et al. (2009) have carried out 
advanced two-dimensional numerical simulations for part of the disk containing the equatorial 
plane by adopting a stratified disk having three layers; cool layer (T ~ 10^ K), warm layer 
(T ~ 10^ K), and hot layer (T ~ 10^ K). Their results show that the magnetic loops have a 
wavelength of ~400 pc and a height of ~350 pc. The typical time scale of the flotation given 
by the ratio H/Va is 10 Myrs for Va of 24 km s~^ corresponding to density of 1000 cm~^ and B 
of 150 /iG (Paper I). These values are fairly consistent with the observed parameters of loops 
1 and 2; the length of the loops, 600 - 800 pc (several times if), and the heights, 200 - 300 
pc (a few times H). In addition, it is important to note that their calculations produced the 
top-heavy loops which show the higher density layer at the top of the loops. This is because 
the loop in the cool layer has initially a small curvature, making the loop top flat, and the 
flat-top loops continuously rise to efficiently accumulate the gas around the loop top. Loop 1 
actually shows that the higher density CO gas covers the lower density Hi is inside the CO 
loop (see Figure 10), indicating that the observations are well reproduced by the local numerical 
simulations by Takahashi et al. (2009). 

Alternatively, the ambient H I gas may be being converted into H2 during the flotation. 
The shock compression by the rising loop in the upper part may lead to form molecular gas 
along the loop which is initially atomic. It is likely that the loop is surrounded by H I gas, and 
the motion of the flotation at Va inevitably causes shocks in the front side. This can lead to 
form H2 in the shock compressed layer along the loop. The time scale of H2 formation is given 
as 10^/n(if) yrs, where n{H) is number density of Hi gas (Spitzer 1978), and is estimated to 
be 10 Myrs for n{H) of 100 cm""^. This time scale is consistent with the present ratio H/Va- It 
is also possible that Hi gas is rising continuously from the disk to the loop top as in the solar 
phenomena. The inside of loop 1 is in fact fllled with lower intensity H I gas. Higher resolution 
H I observations may shed light on this possibility. It is shown by the TRACE satellite at UV 
that such continuous rising loops are actually observed in the sun.^ These actions can also 
accumulate molecular mass in the loop. Takahashi et al. (2009) shows indeed that the flat top 
of the loops continuously rise and accumulates the gas efficiently. The higher density of the 
loop top derived by numerical results support this scenario. 

4-6. Magnetic field 

There is not yet direct measurements by the Zeeman effects of the magnetic field in the 
molecular clouds in the Galactic center region. OH measurements give only upper limits of 
several mC (Killeen, Lo & Crutcher 1992). Observationally, it is a difficult task because the 
splitting is small and the OH spectra is so broad 100 km s~^. 

Indirect arguments on the arched filaments suggest that the field strengths may be as 

^ http://tracc.lmsal.com/ 
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large as 1 mG (Yusef-Zadeh & Morris 1987) and the double helical filaments may also suggest 
the strong magnetic field (Morris, Uchida & Do 2006). On the other hand, the strong field 
may not be occupying the whole volume as suggested by the low-frequency measurements of 
the field (La Rosa et al. 2006). The strong field from 100 fiG to 1 mG is consistent with the 
strong stellar gravity of the central 1 kpc. The strong gravity results in differential rotation 
of the gas with frozen-in magnetic field and strong shearing motion is induced. This shearing 
motion amplifies the field. Thus, the strong field is a natural result of the strong stellar gravity. 

4-7. Comparisons with the Solar phenomena 

The present study revealed helical distributions of the loops. On the solar surface we 
find similar helical distributions which may be explained by MHD instabilities. Figures 7 and 8 
reveal that loops 1 and 2 have helical distributions which change within a small velocity interval 
of ~ 4 km s~^. We shall compare these distributions with solar loops and theoretical calculations 
of magneto-hydrodynamics and discuss formation mechanisms of the helical distributions. 

First, we tried to estimate the amplitude and wavelengths of the helix by applying 
sinusoidal and elliptical functions as follows; 



where Ic and are the center position of an ellipse, a and b are the half lengths of two axes 
of ellipse, A is amplitude of wave and k is the number of waves in the ellipse, O is defined as 
B = arctan {Y{1 — Iq)/ X{b — bo)), and A6 is a offset of G. As a result we find the parameters 
in Table 4 show nice fits to the observations as shown in Figures 7 and 8. The wavelengths of 
helix in loops 1 and 2 are estimated approximately 0.28°(41 pc) - 0.42°(62 pc) and 0.28°(41 
pc), and the amplitudes are estimated approximately 0.2°(29 pc) - 0.34°(50 pc) and 0.26°(38 
pc), respectively. 

Next we compare with the solar phenomena. It has been known that solar prominences 
show helical distributions (Matsumoto et al. 1998). Formation of prominences are discussed 
under the two ideas below; 1) Magnetic fiux tubes having helical shapes originally rise up. 2) 
Magnetic reconnection in the corona forms a helical shape. 

Recent Hinode observations show evidence for 1 (Okamoto et al. 2008) and we discuss 
along the idea 1 and compare the observations with numerical simulations. Matsumoto et al. 
(1998) show that distorted magnetic field lines becoming unstable against the kink instability 
form a helical shape while they rise. In fact soft X-ray observations of solar corona shows such a 
feature; Wu, Huang & Tang (2005) discussed if the helical distributions satisfy the condition for 
kink instability that the number of turns is large enough to cause the instability. The number 
of turns in the helical distrbutions in the top of the loops $ = LioopBg/rBz, where Lioop is the 
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length of the loop, r is the radius of the cross-section, Be and are the axial and azimuthal 
magnetic fields, respectively, are estimated as follows; 
loop 1: $ ~ Stt 
loop 2: $ ~ Svr 

The both satisfies the condition for kink instability $ > 2.5tt (Hood & Priest 1981), suggesting 
that the helical loops are subject to the kink instability of m=l. One possibility of formation 
of helical distributions is the gas in a rotating disk accumulated into a small volume by the 
Parker instability stores the magnetic distortions (Shibata & Matsumoto 1991) and rise up to 
form the helical shape. The present findings of the helical distributions in loops 1 and 2 offers 
a support for that the magnetic floatation is a viable idea to explain these features and more 
detailed numerical simulations are desirable along this line. 

5. Conclusions 

We have presented a detailed study of magnetically floated loops 1 and 2 in the Galactic 
center based on the NANTEN ^^CO and ^^CO(J=1-0) dataset. The main conclusions are 
summarized below; 

1) Molecular loops 1 and 2 have lengths of ~300 - 400 pc and heights of ~200 - 300 
pc projected on the sky with a width of ~30 pc. They are apparently connected with each 
other and have foot points of enhanced intensities and linewidths toward their both ends. The 
size and velocity span of the foot points are 20 - 30 pc and 50 - 100 km s~^, respectively. 
Theoretical works on magnetically flotation loops predict formation of such a loop with the 
foot points as a result of falling-down motion to the disk along a loop and are consistent with 
the observational properties of the loops. 

2) We show that the tops of loops 1 and 2 show helical distributions which may be due to 
magnetic instability. We also note that the current upper limits on the magnetic field strengths 
obtained in OH Zeeman measurements are consistent with the field strength required in the 
present picture ~150 /xG, while direct measurements of the Zeeman splitting yet remains to be 
achieved. 

3) By assuming that the loops are of the same size at a radius R from the center, we 
estimate the geometrical and kinematical properties of the loops such as the projected length 
of a loop, L, the angle subtended by a loop, 6*, and the rotation and expansion velocities, Kot 
and Vexp- The results are as follow; R = 680 pc, L = 350 pc, 9 = 30° and V^ot = 47 km s~^ 
and Vexp = 140 km s~^. 

4) The two loops are associated with Hi gas, which basically shows similar distribution 
and kinematics with molecular gas. H I gas of lower intensity is distributed inside the loops 
where CO emission is not detected. We also identified that the dust emission at 60 fim and 
100 /im obtained with the IRAS is also associated with loop 1. The total mass of the loops 
are estimated by combining the molecular data. Hi data and the dust emission. The total 
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mass of loops 1 and 2 including all these components are estimated to be 2.3(±0.9) x 10^ Mq 
and 3.7(±1.5) x 10^ Mq. The mass of the molecular gas is much larger than the Hi mass, 
indicating that the atomic gas is a major component of the loops. The total kinetic energy of 
loops 1 and 2 are estimated to be ~10^^ ergs for the observed velocity dispersions. 

5) We have examined an alternative idea that supershells are responsible for the loop 
formation instead of the magnetic floatation. We find that the velocity spans of the loops, 50 - 
100 km s~^, are significantly larger than the typical velocity dispersion of molecular supershells 
created by stellar explosions, 20 km s~^, and that there are no known stellar objects inside 
the present loops. The velocity distribution which is dominated by a large monotonic velocity 
gradient is either not consistent with the kinematics expected for a spherical expansion. We 
thus conclude that supershells are not likely cause of the loops. 

6) The recent global numerical simulations of magneto-hydrodynamics in the nuclear disk 
of 2-kpc radius by Machida et al. (2009) indicate formation of the loops is a natural consequence 
of a differentially rotating magnetized gas disk under the strong stellar gravity. By adopting an 
axially symmetric gravitational potential, these authors find that the gas distribution tends to 
obey a m=l mode where half of the disk has lower densities and that the well developed loops 
tend to be formed in this lower density half. In addition, two-dimensional local simulations 
are made by Takahashi et al. (2009), which assume stratified gas consisting of a cool layer 
(T ~ 10'^ K), warm layer (T ~ 10^ K), and hot layer (T ~ 10^ K) to confirm the observations. 
The results show that the flat top loops effectively accumulate the gas and make a dense top 
layer. We note that the properties of the loops shown by these simulations are consistent with 
the observations. We also discuss that the bar-like potential employed to explain the radial 
motions of the molecular gas is compatible with the magnetic flotation. Such a potential may 
be adopted as a next step in the global numerical simulations as a more realistic model. 
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Table 1: Observation summary 
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Tabic 2: Properties of the loops 
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Table 3: Geometrical and kinematical parameters of loops 1 and 2 



Name 


R 
(pc) 


L 
(pc) 


^0 

n 


e 

n 


(km s^-'^) 


(km s^-') 


loops 1 & 2 


670 


360 


27.5 


31.0 


47 


141 


loop3 


1000 


750 


129.5 


43 


80 


130 



Table 4: Fitted parameters of helixes in loops 1 and 2 
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t loop 1 shows two solutions. 
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Fig. 1: Integrated intensity distributions of ^^CO(J=1-0) obtained by NANTEN. Dotted boxes 
in the figures show the regions that loops are located, (a) Velocity integrated intensity distri- 
butions. Integration range in velocity is from -300 to 300 km s~^. Contours are plotted at 10, 
20, 60, 100, 220, 340 K km s'^ in black and 460, 820, 1180, 1540, 1900, 2260, 2620 K km s'^ 
in white, (b) Longitude - velocity diagram. Velocity resolutions are smoothed to 3 km s~^. 
Integrated range in galactic latitude is from -2° to 2°. Contours are plotted every 6 K degree 
from 0.25 K degree (black) and then every 1.2 K degree from 10 the one. 
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Fig. 2: (a, b) Integrated intensity distributions of loops 1 and 2 in the Galactic center. Dotted 
boxes show the region that loops are located, (a) Velocity integrated intensity distributions. 
Integration range in velocity is from -300 to -20 km s~^. Contours are plotted at 10, 20, 60, 
100, 220, 340 K km s'^ in black and 460, 820, 1180, 1540, 1900 K km s'^ in white, (b) 
Velocity - latitude diagram. Velocity resolutions are smoothed to 3 km s~^. Integration range 
in galactic longitude is from -10° to 10°. Contours are plotted every 1.6 K degree from 0.8 K 
degree (black) and then every 3.2 K degree from 6th one (white), (c, d) Integrated intensity 
distributions of loop 3 in the Galactic center. Dotted boxes show the region that loop is located, 
(c) Velocity integrated intensity distributions. Integration range in velocity is from 20 to 300 
km s~^. Contours are plotted at the same levels of figure (a), (d) Velocity - latitude diagram. 
Velocity resolutions are smoothed to 3 km s~^. Integration range in galactic longitude is from 
-10° to 10°. Contours are plotted at the same levels of figure (b). 
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Fig. 3: (a-d) Integrated intensity maps of loop 1 in ^^CO(J=1-0) and ^•^C0( J=l-0). (a) 
Position-position diagram of loop 1 in ^^CO(J=1-0). loop 1 is seen toward / ~ 356° — 358° , 
~ 0° — 1.5° . The integration range in velocity is from -180 to -90 km s~^ . Contours are plotted 
every 10 K km s~^ (black line) from the lowest one of 7.0 K km s~^ {3a) and every 20 K km s~^ 
(white line) from the 6th one. (b) Position-position diagram in ^^CO(J=1-0). The integration 
range is the same in figure (a). Contours are plotted every 4 K km s~^ from 3.5 K km s~^. (c) 
Velocity-galactic longitude diagram of loop 1. The velocity is smoothed at a resolution of 2.0 
km s~^. The integration range in the galactic latitude is shown by the solid lines in figure (a). 
Contours are ploted every 0.12 K deg from the lowest one of 0.20 K degree. An approximate 
position of the 3kpc arm is indicated by a dotted line, (d) Velocity-galactic longitude diagram 
of loop 2. The velocity is smoothed at a resolution of 2.0 km s~^. The integration range is the 
same as figure (c). Contours are plotted every 0.05 K degree from 0.035 K degree. 
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Fig. 3: (e-h) Integrated intensity maps of loop 2 in ^^00(7=1-0) and ^^CO(J=1-0). (e) 
Position-position diagram of loop 1 in ^^CO(J=1-0). loop 2 is seen toward / ~ 355° — 356° , 
6 ~ 1° — 2°. The integration range in velocity is from -90 to -40 km s~^ . Contours are plotted 
every 10 K km s~^ (black line) from the lowest one of 7.0 K km s~^ {3a) and every 20 K km s~^ 
(white line) from the 6th one. (f) Position-position diagram in ^^CO(J=1-0). The integration 
range is the same in figure (e). Contours are plotted every 4 K km s~^ from 4 K km s~^. (g) 
Velocity-galactic longitude diagram of loop 2. The velocity is smoothed at a resolution of 2.0 
km s~^. The integration range in the galactic latitude is from 0.7° to 2.5°. Contours are ploted 
every 0.12 K deg from the lowest one of 0.20 K degree, (h) Velocity-galactic longitude diagram 
of loop 2. The velocity is smoothed at a resolution of 2.0 km s~^. The integration range is the 
same as figure (g). Contours are plotted every 0.05 K degree from 0.035 K degree. 
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Fig. 4: Spectra of loops 1 and 2 in ^^00(7=1-0) with black line and ^^CO(J=1-0) with blue 
line. 
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Fig. 5: The distributions of the peak ^^CO(J= 1 0) velocity in loops 1 and 2. The peak 

velocity is determined by gaussian fitting to the averaged five spectra taken as a cross for each 
position, loops 1 and 2 are ploted 5.2 K km s~^(Bold contours) and 45.2 K km s~^ (dashed 
line). Thin contours are ploted every 10 km s~^ from -200 to -40 km s~^(a), from -120 to km 
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Fig. 6: Velocity channel maps of loops 1 and 2 in CO(J=1-0) integrated over successive 10 
km/s. Contours are plotted every 4.8 K km/s (black contours) from 3.6 K km/s and every 9.6 
K km/s (white contours) from 32.4 K km/s 
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Fig. 6: continued. 
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Fig. 7: Velocity-integrated intensity map helical feature in loop 1. (a) Velocity integrated 
intensity map in ^^CO(J=1-0). The integration range is from -129 to -107 km/s. Contours are 
plotted every 2 K km/s from 2.6 K km/s. Black box indicates the region shown in figure (b). 
(b) Velocity channel maps of helical feature integrated over 4 km/s with interval of 2 km/s. 
Contours are plotted every 0.8 K km/s from 1.6 K km/s. Dotted lines show the fitting result 
with the parameters shown in Table 4. 
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Fig. 8: Velocity-integrated intensity map helical feature in loop 2. (a) Velocity integrated 
intensity map in ^^CO(J=1-0). The integration range is from -86 to -66 km/s. Contours are 
plotted every 2 K km/s from 2.6 K km/s. Black box indicates the region shown in figure (b). 
(b) Velocity channel maps of helical feature integrated over 4 km/s with interval of 2 km/s. 
Contours are plotted every 0.8 K km/s from 1.6 K km/s. Dotted lines show the fitting result 
with the parameters shown in Table 4. 
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Fig. 9: Velocity-integrated intensity map helical feature in loop 2. (a) Velocity integrated 
intensity map in ^^CO(J=1-0). The integration range is from -44 to -34 km/s. Contours are 
plotted every 2.4 K km/s from 2.6 K km/s. Black box indicates the region shown in figure (b). 
(b) Velocity channel maps of helical feature integrated over 4 km/s with interval of 2 km/s. 
Contours are plotted every 1.2 K km/s from 1.6 K km/s. 
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Fig. 10: Integrated intensity maps of Hi and ^^CO(J=1-0) emission. The color image and 
dotted contours indicate Hi emission. Solid contours indicate ^^CO(J=1-0) emission, (a) 
Velocity- integrated intensity maps of loop 1. Integration velocity range is from -180 to -90 
km/s. Hi: Contours are plotted at 42 K km/s. ^^CO(J=1-0): Contours are plotted every 4.8 
K km/s from 3.6 K km/s and every 9.6 K km/s from 6th one. (b) Velocity-integrated intensity 
map of loop 2. Integration velocity range is from -90 to -40 km/s. Contours are plotted at 
the same levels in figure (a), (c) Galactic longitude- velocity diagram of loop 1. The integrated 
galactic latitude range is shown in figure (a) by solid lines. Hi: Contours are plotted at 0.42 
K degree. ^^CO(J=1-0): Contours are plotted every 0.2 K degree from 0.12 K degree, (d) 
Galactic longitude- velocity diagram of loop 2. Contours are plotted at the same levels of figure 
(c). 
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Fig. 11: (a) Velocity-galactic latitude diagram of loops 1 and 2 in ^^CO(J=1-0) integrated 
from 355° to 358° in galactic longitude. Contours are plotted every 0.2 K degree from 0.24 
K (black) degree and every 0.4 K degree from 9th one (white), (b) Velocity-galactic latitude 
diagram of loops 1 and 2 of Hi and CO integrated from 355° to 358° in galactic longitude. The 
color image and dotted contours indicate Hi emission. Solid contours indicate ^^CO(J=1-0) 
emission. Hi: Contous are plotted at 0.42 K degree. ^^CO(J=1-0): Contours are plotted at 
the same levels of figure (a). 
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Fig. 12: CO Integrated intensity distributions superposed on IRAS 100/im. Integration range 
of CO is from -300 to 300 km s~^. Contours are plotted every 15 K km/s from 7 K km/s and 
every 30 K km/s from 5th one. 
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Fig. 13: CO Integrated intensity distributions superposed on IRAS 60/im (left) and 100 /zm 
(right). Integration range of CO is from -300 to 300 km s~^. Contours are plotted every every 
300 K km s"^ from 100 K km s^^ 



34 



ira 3oa '.r-.r Jto ^ Li-si 




ji^tj j^ii :iST List^ Lise :ist jse 



Fig. 14: CO Integrated intensity distributions superposed on 100 fim. (a) Loop 1: Integration 
range is from -180 to -90 km s~^. (b) Loop 2: Integration range is from -90 to -40 km s~^. (c) 
Disk and local components: Integration range is from -40 to 30 km s~^. (d) Loop 3: Integration 
range is from 30 to 200 km s~^. Contours in all figures are plotted every 15 K km/s from 7 K 
km/s and every 30 K km/s from 5th one. 
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Fig. 15: (a) Spatial distributions of column density of atomic hydrogen, N(H), estimated from 
CO and Hi integrated intensity, N(H)gas, vs. N(H) estimated from IRAS 60 fim and 100 /im, 
N(H)dust; assuming dust-to-gas ratio of 100. Black boxes are the plotted regions in figure 8 in 
red. 
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Fig. 16: Scatter plot of the column density of atomic hydrogen, N(H), estimated from CO and 
Hi integrated intensity, N(H)gas, vs. N(H) estimated from IRAS 60 /im and 100 fim, N(H)dust, 
assuming dust-to-gas ratio of 100. Plotted region is the same shown in figure 7. Red points 
show the loop top regions shown in figure 7. Filled triangles denote the foot point of loops 1 
and 2, and filled triangles and diamonds denote the top regions of loops 1 and 3 and top of 
loop 2, respectively. Dotted line show the fitting result for all points; logio(N(H)dust) = 0.76 x 
logio(N(H)gas) + 5.47 
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Fig. 17: Scatter plot of the integrated intensity, W(CO), vs. the column density of molecular 
hydrogen, N(H2), estimated from ^^CO. Solid line show the fitting result; N(H2) = 1.6 x 10^° 
X WC^CO). 




Fig. 18: Schematic image of the position informations of loops 1, 2 and 3 with face on view. 
The parameters in the figure are also summarized in Table 3. 
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